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ABSTRACT 

as : 
as • 

We consider the radiative acceleration to relativistic bulk velocities of a cold, optically 
thin plasma which is exposed to an external source of 7-rays. The flow is driven by 
radiative momentum input to the gas, the accelerating force being due to Compton 
^ | ■ scattering in the relativistic Klein-Nishina limit. The bulk Lorentz factor of the plasma, 

r, derived as a function of distance from the radiating source, is compared with the 
corresponding result in the Thomson limit. Depending on the geometry and spectrum 
{Sj . of the radiation field, we find that particles are accelerated to the asymptotic Lorentz 

factor at infinity much more rapidly in the relativistic regime; and the radiation drag 
is reduced as blueshifted, aberrated photons experience a decreased relativistic cross 
section and scatter preferentially in the forward direction. The random energy imparted 
' to the plasma by 7-rays can be converted into bulk motion if the hot particles execute 

many Larmor orbits before cooling. This 'Compton afterburn' may be a supplementary 
source of momentum if energetic leptons are injected by pair creation, but can be 
neglected in the case of pure Klein-Nishina scattering. Compton drag by side-scattered 
radiation is shown to be more important in limiting the bulk Lorentz factor than the 
finite inertia of the accelerating medium. The processes discussed here may be relevant 
to a variety of astrophysical situations where luminous compact sources of hard X- 
and 7-ray photons are observed, including active galactic nuclei, galactic black hole 
candidates, and gamma-ray bursts. 



Subject headings: gamma-rays: bursts - theory - radiation mechanisms 



1. Introduction 



The ejection of particles by radiation pressure has been considered many times as a possible 
mechanism for producing relativistic outflows from very luminous radiation sources, such as active 
galactic nuclei (AGNs) or galactic compact objects. In the case of an optically thin plasma, the 

institute of Astronomy, Madingley Road, Cambridge CB3 OHA, UK. 

2 Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106-4030. 
3 Department of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599. 



- 2 - 



particle dynamics is dominated by the intense photon field and a flow results from the work done 
by the radiation force on the material. 

In most previous studies it was assumed that the particles scatter photons with a cross section 
independent of frequency and front-to-back symmetric along the incident photon direction, such 
as true in the Thomson limit (Noerdlinger 1974; O'Dell 1981; Phinney 1982; Kovner 1984). Direct 
momentum input to the fluid at rest results when outward-directed photons are removed from 
the anisotropic radiation field, and then scattered away. No net momentum is carried off by the 
scattered photons and the material is accelerated radially. A particle moving outward near the 
source, however, encounters many photons nearly at right angles to its motion, more so as v — > c 
because of aberration effects, and suffers collisional drag by them. Only mildly relativistic flows 
can be produced by an Eddington-limited extended source, the radiation drag actually preventing 
acceleration to very high particle energies (Noerdlinger 1974). In the context of AGNs, the role 
of Compton drag in the dynamics of relativistic jets has been the subject of many studies (e.g. 
Abramowicz & Piran 1980; Sikora k Wilson 1981; Phinney 1987; Melia & K6nigl 1989; Sikora et 
al. 1996). 

In the past few years, the Compton Gamma-Ray Observatory mission has revealed the existence 
of some classes of objects, such as the EGRET (Fichtel et al. 1994), OSSE (McNaron-Brown et 
al. 1995), and COMPTEL (Bloemen et al. 1995) blazars, together with some 7-ray pulsars (Ulmer 
1994), where the hard X-ray and 7-ray fluxes are a significant fraction or completely dominate 
the overall radiation energy budget. This is similar to what sometimes observed by the SIGMA 
experiment onboard of Granat in the case of galactic black hole candidates (Mandrou et al. 1994). 
The large compactnesses inferred from these data naturally lead to theoretical models in which 
Compton scattering in the Klein-Nishina (KN) limit (as opposed to Thomson scattering), together 
with photon-photon pair production may play a role in determining the thermal state and dynamics 
of the source. Such process may also be relevant for studies of the effects of gamma-ray bursts on 
their gaseous environment (their "afterglows"). 

In this paper we study the dynamics of a cold plasma embedded in the "hard" photon field of an 
external source. We adopt a test-particle approach, which is valid only insofar gas pressure gradients 
can be neglected compared to the radiative force and the material is optically thin. We consider the 
model problem of a fully ionized plasma where gravity is overwhelmed by radiation pressure and the 
distribution of particle momenta remains one-dimensional. We ignore stimulated scattering, photon 
absorption, and photon-electron pair production. Free-free absorption is unimportant compared 
to Thomson scattering if <jff = 8 x lCP 47 n e (hv/10~ 3 m e c 2 )~ 7 ^ 2 cm 2 < o~t {hv 3> kT), where n e 
is the electron density. Photo-electron pair production, je~ — > e + e~e~ (threshold hv' > Am e c 2 , 
where hv' is the photon energy in the electron rest-frame) is a higher order process in the fine- 
structure constant, important only at large energies. (The pair production cross section is ~ 30% 
of the KN scattering cross section at hv' = 100 m e c 2 , e.g. Svensson 1982.) In an accompanying 
paper (Thompson & Madau 1999, hereafter Paper II) we show, however, how pair creation induced 
by collisions between hard photons and soft side-scattered photon may increase the acceleration 
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rate by decreasing the inertia per particle, and by increasing the mean momentum deposited per 
scattering. 

The plan of this paper is as follows. In § 2 we discuss the basic theory of radiative acceleration. 
In § 3 we present sample solutions for extended and impulsive radiation sources in the Thomson 
regime. Compton drag by side-scattered radiation is shown to be more effective in limiting the 
bulk motion than is the finite inertia of the accelerating medium. The random energy imparted to 
the plasma by 7-rays can be converted into bulk motion if the hot particles execute many Larmor 
orbits before cooling. The equation of motion in the KN limit is derived and numerically integrated 
in § 4 for a monoenergetic spectrum, under the assumption that the outflowing particles maintain a 
one-dimensional distribution of momenta, i.e. that the plasma remains 'cold' even in the presence 
of recoil. Blumenthal (1974) used the KN cross section to calculate the mean force due to Compton 
scattering on electrons with arbitrary velocity. An approximate solution to the problem of the 
radiative deceleration of a relativistic jet in the KN regime has been recently given by Luo & 
Protheroe (1999). 



2. General formalism 

Consider an optically thin plasma exposed to a photon source. The incident and scattered 
photon momenta (in units of m e c) are denoted by xk and x s k s in the (unprimed) lab frame. Let 
Zr(k) be the specific intensity of the incident radiation in the direction k, J(k) = / I x (k.)dx. In the 
simplest case of a leptonic plasma composed entirely of electrons and positrons, the bulk momentum 
is obtained directly from the normalized phase-space density /(p), 



(173) = / -^/(p)rfV (1) 

J m e c 

Here, T = (1— /? 2 ) -1 / 2 as usual. This equation is easily generalized to include a hadronic component. 
The mean rate of momentum transfer per particle is then given by 



Xsks Xk (l - (3 ■ k)I x (k)dxdadn 



f(p)d 3 p. (2) 



Here, fi is the mean mass per scattering charge, the angular integral is over the solid angle subtended 
by the source, and da is the invariant differential cross section. 

When the number of scattering charges is dominated by pairs, we have assumed implicitly a 
magnetic field strong enough to couple the pairs with the hadronic component of the plasma. The 



minimum flux density needed to enforce that coupling is estimated in §3.5| . One must then describe 
the accelerating particles as a fluid, with 



(3) 
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In a cold plasma composed of hydrogen (of density n p ) and e ± pairs (of density n e +), 

H^m e + m p (l + 2— ) . (4) 
V n P J 

This expression neglects the inertia of the neutralizing electron component, but not of the pairs. 

The Klein-Nishina cross section for scattering of unpolarized radiation is most simply expressed 
in the (primed) electron rest-frame as 



3 fx[ 



i \ 2 r ^' 



x . x s -21 

7 + ZJ " sm X 



_x' s x 

e.g. Rybicki & Lightman (1979). The energy of the scattered photon is 



sin x'd,x'd4>, (5) 



x' 



X ' s l + x^l-cosx') 5 (6) 
in this notation, the scattering takes place through a polar angle x' anci an azimuthal angle <j>' . 

Measuring the angles 9 and 9 S of the incident and scattered photon with respect to the electron 
velocity /3, we can write the following kinematic relations 

x s = x' s T{l + (3 cos 9' s ) 
x = x'T(l + /3 cos 9') 

(1 — p cos 9) 
cos 9' s = cos 9' cos x' + sin 9' sin x' cos 4>' . 

Since the scattering cross section is independent of the azimuth (j)', we must integrate equation 
(j2|) over two angles, x' an d 9'. The evaluation of the integral in equation (Q) will be made by 
assuming that the radiation field is symmetric around the r-axis, and the electrons (and positrons) 
will maintain a one-dimensional distribution of momenta along r, f(pn) = f(p)5(n — f). Note 
that multiple scatterings can be safely neglected, as the singly scattered photons are beamed into 
a cone of half-angle ~ 1/T along the electron direction of motion. Further scatterings are therefore 
suppressed by the (1 — (3 ■ k) ~ (2r 2 ) -1 factor in the momentum transfer equation. Single side- 
scattered photons will, however, be absorbed via 7 + 7 — > e + e~ if the compactness of the 7-ray 
source is high at x > 1. We defer a detailed discussion of the effects of pair creation to Paper II. 



3. Thomson scattering 

To provide some physical insight and a framework to interpret our numerical results in the 
relativistic limit, we shall first discuss the acceleration of a test charge in the Thomson scattering 
approximation, x' <C 1. In this regime, x' s ~ x' by equation (||). Thus, in the electron rest-frame, 
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photons of unchanged energies are scattered into an angular distribution that is forward-backward 
symmetric along the incident photon direction. 

A particle accelerated by an external photon source will, in general, follow a curved trajectory. 
The problem simplifies when the source is cylindrically symmetric. From equation (|2|), a parcel of 
matter moving along the axis of symmetry with speed j3 = c~ 1 dr/dt is accelerated at the rate 



d{T(5) _ a T r 2 



dt fjLC 



2 



{l-(3f / 1(6) cos 8dQ- (3 / 7(0) (1 - cos 8) 2 dfi 



(8) 



The first term in brackets represents the collimated flux of redshifted photons which accelerate 
particles radially, while the second term represents the dragging effect of the isotropic component 
of the radiation field.[] The right-hand side of equation (|J) vanishes for T = T eq . Particles are 
accelerated away from the source as long as V < T eq . If T > r eq , the force reverses, being now 
directed inward. This saturation of the Lorentz factor is due to the aberration into the forward 
hemisphere of blueshifted photons, as seen in the electron rest-frame. A particle will follow the 
r = T eq equilibrium trajectory (zero-inertia limit) until, as it gets farther from the source, its 
acceleration is effectively limited by the large Doppler shift to the red and by the r~ 2 decline of 
the radiation flux. 



3.1. Uniformly radiating sphere 

Consider a sphere of radius R and uniform brightness 7. Close to the extended source, the 
equilibrium Lorentz factor is 

(9) 

no matter how high the intensity 7, or whether the inertia is dominated by hadrons or leptons. In 
the point-source limit (R — > 0) the photons stream radially. For a steady flow, equation (||) reduces 
to 

d ^-=T\i-pr4. (io) 



The parameter t is the dimensionless compactness, rescaled by the inertia per scattering charge, 
i = £(m e /fi) = La T /(^fj,c s R). It may be expressed as i = 0.5(m p / n)(L/ Le)(Rs / R) for a source 
of mass M, where Le = 4TrcGMm p /aT is the Eddington luminosity and i?g = 2GM/c 2 the 
Schwarzschild radius. For a motion starting at radius r, the Lorentz factor attained at infinity can 
be obtained by solving the algebraic equation 

2r3 s (i + /? 3 )-3r ps + i = 3^, (ii) 



4 This term vanishes in the case of a point-like source, /(k) = F(r)S(k. — r), where F(r) is the radiative flux at 
radius r. 
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which, in the limit T ps > 1, yields T ps ~ [3£R/(4r)]^ 3 = [3Lm p R s / '(SL^r)] 1 / 3 . Effective accel- 
eration of an e~p plasma requires a compactness that is larger by a factor m p /m e than for a pure 
pair plasma, as a conseguence of the greater inertia per unit cross section. 

Matching this solution with expression (^) yields the asymptotic bulk Lorentz factor F 

£1/4 

(Noerdlinger 1974). The cross-over radius beyond which the non-radial component of the radiation 
field causes negligible drag is r c ~ 0.8^ 1//4 i?. Larger values of can only be obtained if the particles 
are injected with relativistic bulk velocities at r > r c . Bulk motion starting with T(r) ^> r eq (r) at 
a distance r < r c from the source is quickly decelerated to V ~ r eq , and the excess kinetic energy 
converted into a collimated beam of upscattered photons. 



3.2. Non-uniformly radiating disk 

Infinite, Keplerian accretion disks around black holes provide an example of non-uniform ra- 
diating sources where Compton drag is most severe. The surface flux distribution as a function of 
equatorial radius R is given by (Shakura and Sunyaev 1973) 

F(R) = [l - {m s /R)^] , (12) 

where M is now the black hole mass, and M is the accretion rate. The peak flux, -F max = 3GMM / 
(56-71-/?!^), is reached at i? max = (49/12)i?s. One can understand the essential features of the 
accretion disk model by approximating the flux distribution as 

F(R) = \o if E<iw (13) 

In this case the equilibrium Lorentz factor increases only 

increase for a spherical source). In the point-source limit, one derives Tp S ~ {0.751 R^/r) 1 / 3 , where 
£ is now the nominal disk compactness, £ = 3GMM&T / (28-71 /ic 3 'R^ ax ) . Again, by matching the 
two limiting solutions, the asymptotic bulk Lorentz factor is found to scale as ~ £ ' (Kovner 
1984; Phinney 1987). The exact solutions for a spherical uniform source and a Keplerian accretion 
disk are shown in Figure 1. The latter is a much less efficient accelerator than the former because 
of the large photon drag associated with a radiation source extending to infinity. 



3.3. Compression of the accelerating medium 

The radiative force acting on an external medium induces a flow that is strongly time-dependent. 
One important effect of acceleration to bulk Lorentz factor T is the bunching of the accelerated 
material into a shell of thickness Ar ~ r/T 2 . Let us now solve for the density of the accelerated 
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material. The equation of continuity of the hadronic component is 



Material initially at radius r$ reaches a radius 

r l ff 3 d,B' 

r = r + c P(t')dt' = r + c P' -f- (15) 
Jr /c JO J\p') 

at time t (speed /?). Here, dp/dt = f(P). The total duration of the acceleration is 



ro 



P dp' 



c Jo f{P') 

Eliminating ro from these two equations and taking d/dr at constant time t, one derives dp/dr = 
f(P)/c(l — P), and we have 

1 d r ,. 1 dp 



Tn' p dV p ' l-pdt 



This integrates to 



^ = -!-~2r 2 (r»i), (is) 

n p0 1 - P 

where n p o is the proton density in the undisturbed ambient medium. 

There is a corresponding amplification of the non-radial component of a magnetic field en- 
trained in the flow, 

A = K = ^_. (19) 

Bo n p0 1 - P 



3.4. Thin photon shell 

Gamma-ray sources such as blazars and gamma-ray bursts are impulsive. It is instructive, 
therefore, to consider a photon source that maintains a constant luminosity for a time At <C r/c. 
This source can be visualized as a uniform shell of radius r = ct, thickness cAt, and total flux 
F = J F x dx. We focus here on the effects of a radiation pulse and take the photons to move 
radially before scattering. 

The compactness within the shell can be expressed directly in terms of the flux of 7-rays, 

t(-) ^ra T -^3, (20) 
and thence in terms of a characteristic photon optical depth,[] 

T c = crj^At — — n, (21) 

m e c z 

5 This quantity is in fact the optical depth <JTnj(xbr)cAt multiplied by the frequency Xhi at which the spectrum 
peaks. 
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such that £(R/r) — » r c (r/cAt). Material overtaken by the photon shell at radius r, and accelerated 
from rest, develops a large bulk Lorentz factor when T c (m e //i) 3> 1. Integrating equation (|K]) in 
this limit yields 

3 3jRAr 3 /m e \ / Ar \ 

at a radius r + Ar. As the difference between the particle speed and the speed of light decreases 
as (Ar) _2//3 , the accelerated matter will surf the photon shell over an extended range of radius, 
Ar = §r 2 cAt. The maximum Lorentz factor is therefore 



Tc ( m e 
2 { /i 



(23) 



p<^^x^^-7^r 1 fM7^- 1 gcm -a. (24) 



when the shell is thin enough that the acceleration length is much less than the radius, Ir^^cAi <C 
r (see Figure 2). Only a limited amount of material can be accelerated to this value. Equating the 
kinetic energy within a volume (47r/3)r 3 with the energy of the photon shell implies a very high 
pre-existing optical depth to Compton scattering, tt = OTprjm v ~ 6. Since our calculation of 
Tmax is valid only for small tt, the ambient density at radius r is bounded above by 

3 

cr T r~""^ \300y \10s7 \m p ) v 2r 2 iax cAt / 

Note the strong dependence on /i. We will show is Paper II that, even if the ambient medium 
is composed of a hadronic plasma initially (as is expected in most cosmological gamma-ray burst 
models), the mean mass per scattering charge will be rapidly reduced to \i ~ m e by pair creation. 

Sidescattered photons will provide additional drag that reduces T below the value given in 
equation (^3[). As the photon shell propagates ahead of an accelerating parcel of matter, it generates 
an isotropic radiation field of photon density 

r 2 . s 

- n s a T ~-, (25) 



IR 



where n s = n p + 2n e + is the density of scattering charges in the ambient medium at rest, and 
r 2 /£R is the lengthscale over which isotropic scattering occurs (before the scattering medium gets 
accelerated to relativistic velocities). The net effect is to introduce a second, negative term into 



equation (10), which becomes 



/ 1 4 - r 2 \ IR 
— = — - -T 2 n s a T ^- -^t 26) 
dr \^T 2 3 £R) r 2 y J 

at large T. Acceleration up to Lorentz factor r max is then possible only if the ambient density is 
lower than 

P < TF^l a7- V 2 ') 

!6r^ ax (T T cl\t 

This bound is stronger by a factor ~ T~ l than equation p^): in the Thomson scattering regime, 
the bulk Lorentz factor is actually limited by side- scattered radiation. This conclusion changes when 
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the photon source is very hard, because side-scattered photons are immediately consumed by pair 
creation (Paper II). A compact source of 7-rays always tends, in this sense, to maintain a radial 
photon distribution. 



It is well known that a hot (relativistic) plasma in an anisotropic radiation field feels a much 
larger radiation pressure than a cold (non-relativistic) gas of the same inertial mass (O'Dell 1981), 
and tends to drive itself away from the radiation source with momentum derived largely from 
the anisotropic loss of its own internal energy. This "Compton rocket" is, however, ineffective at 
generating bulk relativistic motion from internal random kinetic energy, as it is always accompanied 
by catastrophic Compton cooling (Phinney 1982). Nonetheless, by integrating the momentum 
equation in the limiting case of a cold plasma we are actually underestimating its bulk Lorentz 
factor: the gas will be continuously heated by the hard radiation field itself, hot particles will 
radiate most of their energy in the direction of the radiation source, and part of the energy of 
relativistic random motion will be converted into bulk motion. We term this supplementary source 
of momentum the 'Compton afterburn'. The afterburn can be largely neglected in the case of 
pure KN scattering without pair creation (§4.1). (In Paper II we will consider a situation in 
which relativistic e + e~ pairs are injected into an anisotropic photon beam by collisions between 
side-scattered photons and the main beam - thereby inducing a net force on the plasma as they 
Compton cool.) Nonetheless, we discuss it here for completeness. 

The afterburn is strongest when the radiative force on the relativistic, cooling particles is cou- 
pled to the rest of the fluid through a background magnetic field. Above a critical flux density which 
we estimate, this allows the colder particles to dominate the inertia of the flow, and the momentum 
gained by Compton scattering to be more effectively absorbed rather than being radiated away. 
The afterburn is easiest to analyze when the period of a cyclotron orbit is much shorter than the 
Compton cooling time. Then the linear momentum imparted to one charge is effectively shared 
with the others before the charge loses its excess energy in the bulk frame. We also assume that 
the acceleration time is long compared with the Compton cooling time, so that the proportion of 
energetic particles is small. This is the case when the energetic charges are continually regenerated, 
e.g., by pair creation. Finally, we assume that the magnetic field is non-radial (as is appropriate to 
the rest frame of a plasma that has been accelerated to relativistic speed). We divide the velocity 
of the charge into components (3 e ±_ and /3 e ii perpendicular and parallel to B, and denote by 77 the 
angle between the cyclotron motion and the radial direction. 

The charge is immersed in a collimated photon beam of energy flux 



3.5. Compton afterburn 





(28) 
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in the boosted frame of the bulk flow. It feels a radial force (O'Dell 1981) 
d ( Pr \ CTtFy 



dt \ m e c) m P c 2 



1 - (3 e ± COS if) 7 e (1 - P e ± cos 1])(-(3 e ± COS 7/) + 1 . (29) 



The net radial momentum imparted to the cooling charge is obtained by averaging over a cyclotron 
orbit, 77 — > 77 + 2-7T, 

At the same time, the net energy lost to Compton cooling is 

cfye \ o- T F r 2 / 1 



(it / TOpC 



7e Q/3e±+/3e)- (31) 



When the energetic charges are created with /3 e u = (e.g., pair creation between a hard radial 

photon and a soft side-scattered photon), the cooling rate is dj e /dt = — (3o"Ti ? r/2m e c 2 ) 7 f 01 and 
the radial momentum gained as the particle cools is 

d(p r /m e c) _ _ 2 7 2 

die - 3( 7 2 -l)- {d2) 

To leading order in 7e , 

/V ' ^- (33) 



m e c 3 

In this case, the source of energy for the charge is external to the accelerating medium. The charge 
carries an excess lab- frame momentum (2T)'j e m e c 1 which is supplemented by (33) to yield a total 

§< 



5 r(2r)7 e m e c (in the limit 7e ^> 1). 



These expressions change slightly when the energetic charges are distributed isotropically, and 
reduce to those found by O'Dell (1981). The rate of Compton cooling is d^ e /dt = —(4:aTFr/3m e c 2 ) 
7 2 /3 2 , as usual, and 



d(p r /m e c) + \ 



die 2( 7 | - 1) 

The net radial momentum acquired is 



(34) 



o7e, (35) 



no 



m e c 2 

only 3/10 of the value obtained for charges that are created moving parallel to the photon beam. 

These calculations assume that a high energy particle will execute many cyclotron orbits before 
Compton cooling. Let us check when this assumption is valid. The rest frame magnetic flux density 
is related to the initial (non-radial) field before acceleration by 



(36) 
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(eq. Each fresh (radially-moving) particle has a rest-frame energy 7 e m e c 2 ~ E/2T, and so it 
takes a time 

, JeTTleC 1 -E 

^ c = ~lB r ~ 4f2 (37) 
to complete one radian of a Larmor orbit. The Compton cooling time is 

j 2ui e c 2 
6^ e aTrr 



The energy flux in a thin photon shell is related to the characteristic optical depth r c (eq. [21J ) 
through aTF/m e c 2 = t c / At, and so the ratio of these two timescales is 



^ ool 32r 5 \m e (?) \eB AtJ T 5 \m e c 2 J \3 x W' 6 G J \10 s 

The compression of the accelerating medium has a related consequence: the mean energy 
(7e) of the scattering charges in the bulk frame increases due to adiabatic heating, at a rate 
(7e) _1 {d(le) I dt) = ^n'~ 1 (dn' /dt) = ^T~ 1 (dF/dt) in the lab frame. A balance between compres- 
sional heating and Compton cooling results in an equilibrium energy (7 e ) that is only mildly rela- 
tivistic. As a result, we will assume that the heating process has isotropized the momenta of the 
pairs in the bulk frame: 

( 7e ) dt 3Tdt 3t° ool ' 1 ] 

The reference (lab-frame) cooling time is i° ool = m e c 2 / cjt^ Ft — (rxb r /2r c )At. The acceleration 
rate is decreased by a factor ~ (7e) _1 from the value calculated above for a cold plasma: 



■'accel 



dT/dt 

^From equation ([U]), we have 



(7e)4cel- (41) 



r accel Zi r 1 

at T 3> 1. Even in the case of a pair- loaded plasma (fi ~ m e ), one sees that energetic particles cool 
faster than the bulk flow accelerates. The same conclusion holds for scattering in the Klein-Nishina 
regime (§4), but not when pair creation drives the acceleration (Paper II). 



4. Klein-Nishina regime: monoenergetic spectrum 

When x' > 1 recoil can no longer be neglected. In the general case, the integral over the 
scattering angle \' °f the Compton force on particles moving along the r-axis, 



d 1 f x q cos 9 q — x cos 8 

— (TP) = ; 

at \ic 



(Tp) = --^ / " 5 ™ (1 - f3 cos 9)I x dxdadn, (43) 
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can be performed analitically by making the variable change x' ~ > x ' s i an d integrating the resulting 
polynomial in x' . The mean rate of momentum transfer can then be written as 



Urn = - T 

at lic 



(TP) = -L I (l - f3 cos e)K(x') 



. T(cos6-l3) 
cos 9 + 



x(l — (3 cos 

where K(x') = J x' s (l — cos x')d(r/(TT can be evaluated in closed form, 



I x dxdn, (44) 



K(x'\ 1 



Ax' 2 



x 



n 



2a/ -3, ,„ „ t\ -10x /4 + 51x' 3 + 93x /2 + 51x' + 9' 
m(l + 2x ) + 



2x' v 1 3(1 + 2x') 3 



(45) 



(cf. Blumenthal 1974). This function has limiting behavior 



K{x') « " 5 /6 + ••••] (46) 

oar 

for x' 1, and 

K(x') w x'(l - 21x75 + ) (47) 

for x' <C 1. The first term in the last expansion, when substituted back into equation (^4|) in the 
limit x' <C 1, reproduces the standard Thomson results of equation (|8|) (for li = m e ). ^] 

Figure 3 shows the integrand function in equation fl44]), P = (1 — (3 cos 6) K (x'){cos 9 + T (cos 8 — 
P)/[x(l—(3 cos^)]}, versus incident photon angle for different photon energies and particle valocities. 
The aberration effect of outward-directed photons into the forward hemisphere (as observed in the 
electron rest-frame) is clearly seen, together with the reduced radiative force as 7-ray photons are 
preferentially scattered in the forward direction. In the Thomson limit, photons at angles cos# < j3 
with respect to the direction of motion are seen in the rest-frame of the electron as blueshifted and 
inward directed, and work to decelerate the flow (P < 0). The same criterion does not apply in 
Klein-Nishina, as blueshifted, aberrated photons experience a decreased relativistic cross section 
and are scattered less efficiently. When (3 = 0.8, for example, the rate of momentum transfer 
becomes negative at angles cos 9 < 0.71,0.35, and 0.13 for photons with energy x = 0.5,3, and 10, 
respectively. 

The complexity of the KN cross section foils analytic calculations of the particle trajectories. 
We have numerically integrated equation ( f44"l) with the Runge-Kutta method assuming a monoen- 
ergetic photon spectrum. The equilibrium Lorentz factor as a function of distance from a radiating 



6 The force on a e~p plasma at rest relative to a point-source of mass M and monoenergetic flux F(r) is, from 
equation ([li]), d(T/3)/dt = (<jt /m p (?)F(r) K(x)(l + x) jx. The Eddington luminosity for which this force just balances 
the gravitational attraction at all radii r is then given by 

x 

K(x)(l + x) ' 



T 4ivcGMm p 
Le = 



where the term in parenthesis (equal to 2.4, 5.9, and 11.8 for x — 0.5, 3, and 10, respectively) is the relativistic 
correction factor to the classical expression in the Thomson limit (Blumenthal 1974). 
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sphere and a disk is shown in Figure 4 for different energies of the incoming photons. While in 
the Thomson limit the radiation drag imposes a linear relation between r eq and distance r from a 
spherical source, in Klein-Nishina the importance of this braking force is reduced: the zero-inertia 
limit is reached then at higher particle velocities. The net momentum transfer in each interaction 
with a 7-ray photon is significantly lower, however, and it is actually more difficult (i.e. larger com- 
pactnesses are needed) for an electron to approach the zero-inertia limit (Figure 4). In Figure 5 the 
asymptotic Lorentz factor at infinity is depicted as a function of source compactness for different 
incoming photon energies. A comparison with the Thomson limit shows that 7-ray photons with 
x = 10 (say) are significantly more efficient at accelerating particles only for compactness £ > 200, 
and more so in the extended disk case (i.e. when the isotropic radiation field is larger) than for a 
spherical source. In the general case, the detailed dynamics of an the plasma will depend on the 
energy spectrum of the incoming radiation (Paper II). 

The KN reduction of the photon drag can be easily quantified in the extreme case of a jet 
moving with relativistic speed towards a point-source of 7-rays. From equation (^4|) we derive 

>>-^+Mff^)K)]- 

where x' = xT(l + 0). The term in parenthesis, equal to 0.094, 0.016, and 0.006 for T = 2 and 
x = 0.5, 3, and 10, respectively, is the relativistic correction to the braking force in scattering 
Thomson. 



4.1. Compton heating 

Together with momentum, energy will be transferred from the radiation field to the gas at a 

P cos 6)I x dxdadil 

I x dxdfi. (49) 



rate 

/dE\ = 1_ f x s -x 



\ dt I Lie 2 J x 

r/3(cosfl-/3) 



^ / (1-3 cos 9)K(x') 

IXC 1 J 



1 + 



x(l- PcosO) 



Equations (44) and ( f49| ) assume that all the particles in a small volume element of gas (that moves 
with bulk Lorentz factor V) are at rest with respect to that volume element. Because of recoil in 



each collision, however, the plasma will be Compton heated (hence (dE/dt) ^ dT /dt in eq. 4£) to 
relativistic temperatures, and in general these equations should be integrated over the appropriate 
distribution function. The integral can can solved numerically for a momentum distribution that 
is isotropic in the electron rest-frame (e.g. a relativistic Maxwellian) provided the relaxation time 
of the plasma is short compared with the dynamical time. 



7 In the absence of any other cooling or heating mechanisms, a thermal gas will be driven towards the Compton 
temperature Tc = m e c 2 (x) /(4fc). For a power-law spectrum with energy indices a — and —1.5 below and above 
the break frequency Xbr = 1 (characteristic of gamma-ray bursts), one derives (a;) = 0.74 and T c = 1.1 x 10 9 K. 
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In the case of a cold plasma at rest in an isotropic photon bath, the Compton heating rate per 
unit volume can be written, from equation ( |49| ) with x' = x, 

f \ K ( X )1 < X 

H = &im s OT I xl x dx, (50) 

J I x _ 

where the term K(x)/x (equal to 0.28, 0.04, and 0.008 for x = 0.5, 3, and 10, respectively) is the KN 
reduction factor to the classical Thomson formula. It has been recently pointed out that Compton 
heating of electrons by hard X-ray background photons may provide a significant energy source for 
the intergalactic medium (Madau & Efstathiou 1999). It is crucial in the cosmological context to 
use equation ( |50| ) rather the Thomson limit, as relativistic corrections become increasingly more 
important at early times when the peak in spectral power of the X-ray background (observed today 
at ~ 30 keV) is blueshifted to higher energies. 

It is also worth noting that the efficiency of the Compton afterburn should be further dimin- 
ished in the KN regime, as the decreased relativistic cross section reduces the anisotropy (between 
highly scattering electrons moving towards the source and those, less scattering, moving away) that 
ultimately drives the effect (Phinney 1982; see also Renaud & Henry 1998). The afterburn may be 
a supplementary source of momentum if the energetic particles are injected by an external energy 
source, and if they execute many Larmor orbits while cooling (§[3^5|). However, in the present con- 
text (in which the particles are assumed cold before scattering) only a small increase in the bulk 
momentum will result from the afterburn effect. Even if we suppose that all the energy transferred 
from the radiation field to the gas is actually converted into bulk motion, and integrate equation 
(|49| ) with (dE/dt) = dT/dt, we get asymptotic Lorentz factors that are at most ten percent higher 
than those computed from the momentum equation. 
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Fig. 1. — Left: Bulk Lorentz factor of a test particle as a function of distance r from a uniformly 
radiating sphere of radius R and compactness t. The equation of motion has been integrated in the 
Thomson regime assuming the particle is initially at rest. From top to bottom: I = 1000, 100, 10 
(solid lines). Dashed line: Zero-inertia limit. Right: Same but for an infinite, Keplerian accretion 
disk (see text for details). 
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Fig. 2. — Bulk Lorentz factor of a test particle (// = m p ) as a function of distance r (in cm) from 
a point-like source of luminosity 10 52 ergs s^ 1 . The equation of motion has been integrated in 
the Thomson regime assuming the particle to be initially at rest. Solid curves: radiation pulse of 
duration 1 sec. The radiative force vanishes when the photon shell moves past the particle. Dashed 
curves: steady source. 
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Fig. 3. — Specific rate of momentum transfer per unit solid angle (see eq. 44), P = (1 — 
j3 cos 9)K{x'){cos 9 + T {cos — (3)/[x(l — (3 cos 9)]}, versus incident photon angle for different photon 
energies x (in units of m e c 2 ) and particle valocities f3 (in units of c). Solid line: x = 0.5. Dashed 
line: x = 3. Dash-dotted line: x = 10. Dotted line: Thomson limit. 



-19- 








i ■ i 



i ■ i_ 



J_ 



' ' ' 



0.2 0.4 0.6 
log r/R m£ 



0.8 1 







~r~r 



0.2 0.4 0.6 
log r/R m£ 



0.8 1 



Fig. 4. — Top left: Equilibrium Lorentz factor (zero-inertia limit) of a test particle as a function of 
distance r from a uniformly radiating sphere of radius R. The photon spectrum is monoenergetic 
and the equation of motion has been integrated in the Klein-Nishina regime. Solid lines: Values 
obtained for different energies of the incoming photons: x = 10, 3, 1 (from top to bottom). Dotted 
line: Thomson limit. Top right: Bulk Lorentz factor as a function of distance. The equation of 
motion has been integrated using the relativistic cross section and assuming the particle to be 
initially at rest. The radiation spectrum is monoenergetic with photon energy x = 10. From top to 
bottom: I = 1000, 100, 10 (solid lines). Dashed line: Zero-inertia limit. Note how, in the relativistic 
limit, the particle velocity saturates to its asymptotic value at infinity much closer to the radiation 
source than in the Thomson regime (cf. Fig. 1). Bottom: Same for a Keplerian accretion disk. 
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Fig. 5. — Asymptotic Lorentz factor at infinity as a function of compactness for a particle initially 
at rest. Solid line: KN cross section, incoming photon energy x = 3. Dashed line: Same for x = 10. 
Dotted line: Thomson limit. Left: Uniformly radiating sphere. Right: Keplerian accretion disk. 



